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On parton distributions in a photon gas 
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Abstract. In some cases it may be useful to know parton distributions in a photon gas. This may be 
relevant, e.g., for the analysis of interactions of high energy cosmic ray particles with the cosmic microwave 
background radiation. The latter can be considered as a gas of photons with an almost perfect blackbody 
spectrum. An approach to finding such parton distributions is described. The survival probability of ultra- 
high energy neutrinos traveling through this radiation is calculated. 
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1 Introduction 

There are reactions in which photons manifest hadronic 
properties [T] . In an analysis of such reactions it may be 
possible to employ formalism used for the case of hadrons, 
at least partly. For example, one can introduce the struc- 
ture function characterizing the parton densities in the 
photon [2,3J. It was suggested in the early 1970's that 
information on this structure function might be accessi- 
ble by deep inelastic electron-photon scattering at e^e~ 
colliders [HIS]. Theoretical investigations at that time al- 
ready revealed the basic logarithmic dependence of the 
function on the four-momentum transfer squared [F', 
[7j. Quantum-chromodynamics (QCD) corrections to the 
pointlike structure of the photon were calculated . Evo- 
lution equations for the parton densities as well as the 
properties of the corresponding solutions were under scrutiny 
(see, e.g., |9]). A description of recent developments and a 
more complete list of references can be found, for example, 
in[ll[lin]. 

The universe is filled with the cosmic microwave back- 
ground (CMB) radiation which can be considered as a 
gas of photons with an almost perfect blackbody spec- 
trum at temperature T « 2.725 K [12]. The CMB plays 
an important role in astrophysics, for example, providing 
a medium which inevitably interacts with high energy cos- 
mic ray particles [l3 l fT4 ] fT5] . We have recently pointed out 
that owing to the existence of this radiation the parton 
content of the real photon may also find a non-trivial as- 
trophysical application 16J. It may be therefore useful to 
know parton distributions associated with the CMB. 

The typical energy of a photon of the CMB at the 
temperature 2.725 K is very low, about 10~^ eV. Never- 
theless one should not be embarrassed by such energies. In 
fact, not every reference frame is suitable for an approach 
based on the parton model. There are frames in which the 
parton densities are undefined, as for example in the rest 



frame of a particle whose partonic structure is studied. 
One can, however, successfully use the parton model in 
a coordinate system with respect to which the particle is 
moving with very high momentum (this question is nicely 
discussed, e.g., in [TT]). Thus, a reference frame fixed to a 
high energy cosmic ray particle propagating through the 
universe and interacting with the CMB may be an appro- 
priate one. 



2 Derivation of parton distribution functions 

Let us consider a photon gas with blackbody spectrum. 
The number of photons in volume V at temperature T 
with energies between a; and iv + dw is given by [IB] 



N{Lu,T)duj 



V 



(1) 



where k is Boltzmann's constant, c is the speed of light 
and h is the reduced Planck constant. Henceforth we set 
c = h = 1. Note that ([1]) is a Lorentz invariant. 

Dividing ([T]) by the total number of photons 



(2) 



(where ((s) is the Riemann zeta- function [19]) yields 
the volume-independent probability to find a photon in 
this energy interval 



n(tj, T)dLL} 



u duj 



2C(3)fc3T3 e^A"^ - 1 ■ 
One can see that ([3]) is only a function of the ratio 



(3) 



UJ 

kT 



(4) 
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and may be rewritten as 
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n{x)dx 



X dx 



2C(3) 



(5) 



In order to make the subsequent discussion more cus- 
tomary, let us regard, just formally, n{x) as the probability 
density function to find a photon carrying the "fraction" x 
of the energy kT. We emphasize that x may vary, in the 
range from zero to infinity (this is the reason of using the 
quotation marks). 

We define the sought-for parton distributions as func- 
tions of this dimensionless variable x in the following way: 



g(x,g2) 



y 



n{y) gi 



y 



(6) 



^{^iQ'^) parametrizes the probability density to 



here q 

find a quark of flavor q with the momentum fraction ^ 
in the photon probed by a hard scattering at virtuality 
scale Q^. 

Though x and y in general vary from zero to infinity, 
their ratio x/y in ([6]) is always less than or equal to unity 
and the function q''{i,, Q^) retains its usual meaning. We 
adopted the result of the quark-parton model [5] : 



8^(1-^)-! 



+ (1-0'] In 



, (7) 



where A''^ is the number of the quark colours, a is the 
fine structure constant, Cq and rrig are the electric charge 
and mass of the quark g, respectively. 

Using (O and we have numerically solved ([6]) for 
u and d quarks. We set Nc = 3, = 10 GeV^ and niu — 
TTid = 0.2 GeV. The obtained results are depicted in Fig.[TJ 
To make the picture more complete, the distribution of the 
photons n{x) is also displayed. One can see that it is more 
probable to find the quarks with small values of x, while 
the number of photons decreases as x approaches zero. 
The corresponding antiquarks behave analogously. 

In principle, more precise calculations may be per- 
formed when needed, for example, by using more real- 
istic parton densities instead of ([7]). The possibility of a 
presence of strange and heavier quarks in the photon can 
additionally be taken into account. But the main task of 
this paper is to give a qualitative description of the ap- 
proach itself. For this reason we did not show the Q^- 
dependence of the parton distributions in detail (this is 
weak, being logarithmic, and does not qualitatively change 
the results). 

We note the Lorentz invariance of this formalism which 
follows from the fact that ([1]) and ([7]) are defined in a 
Lorentz-invariant way. 

One may also derive distributions of the charged lep- 
tons fi{x, Q^) exactly as above, just by making the follow- 
ing changes in ([7]): ^ 1, — > 1, to, — * m; [5], where 
the subscript I refers to the leptons. 




Fig. 1. The parton distributions in the photon gas as functions 
of x: u quark (dashed curve) and d quark (solid curve). Note 
that is fixed at 10 GeV'^. The distribution of the photons 
n{x) is also displayed by the dotted curve 
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Fig. 2. Diagrams illustrating a the inclusive on-shell W bo- 
son production !/e7 — > W^X; b charged current neutrino scat- 
tering off quarks (antiquarks). In this paper we take only the u 
and d quarks (antiquarks) into account and neglect Cabibbo- 
Kobayashi-Maskawa mixing 
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[24] . Following them we will calculate the damping rate of 
the neutrinos in the CMB. 

It is convenient to begin with a consideration of the 
interactions at the subprocess level explicitly writing down 
some kinematical quantities. Thus, let (t(s, Q^) be the to- 
tal cross section of a neutrino-parton scattering, where 
standardly 



s = 2Ee(\ - cos( 



(10) 



here E and e are the energies of the incident neutrino 
and parton, respectively (their masses are neglected), 9 
is the angle between the momenta of the colliding parti- 
cles. Note that is the same at the neutrino-photon and 
neutrino-parton levels [25] . 

For a static target the damping rate 7 would simply 
read 7 — cm (recall that c = 1) with n and a being the 
number of target particles per unit volume and cross sec- 
tion of a given reaction, respectively. In our case, however, 
the scatterers are not at rest and we have to take their en- 
ergy distribution into account. We define the damping rate 
at fixed and T in the form 



Fig. 3. Dependence of the survival probability of an UHE r r d'^o; 2q'^(e/L0 Q^) 

electron neutrino emitted at Zs = 10, 15 and 20 on its present 7(-^) ^ '^^ 7 73 oVfeT 1 — '^{s,Q ), (11) 



day energy E. The absorption channel is v^'y W X ji^^ 



3 Neutrino absorption by the CMB radiation 

Here we illustrate how the model presented in the previ- 
ous section can be applied. For this purpose we consider 
two possible channels of absorption of ultra-high energy 
(UHE) neutrinos by the CMB radiation, namely the in- 
clusive resonant boson production 

Vel W+X (8) 

and the neutrino-photon deep inelastic scattering 

^ e-X. (9) 

The corresponding diagrams are shown in Fig.[2| A chan- 
nel closely related to (Uj) has been proposed in [20] . Let us 
discuss only the case of the electron neutrino since the 
other flavours can be treated likewise. 

It is obvious that the incident neutrino has to possess 
extremely high energy to interact with the CMB photons 
via these channels [16] . This is similar to the situation 
when one needs to probe the region of small values of the 
Bjorken variable in the nucleon. As objects of QCD its 
valence quarks are surrounded by a cloud of soft virtual 
gluons which may generate quark-antiquark pairs domi- 
nant at ^1. The lower values of are to be studied, 
the higher energies in the center-of-mass system of the col- 
liding nucleons are required. Analogously we regard the 
CMB radiation as such a cloud whose quark component 
can be resolved by the UHE neutrinos. 

There have been works devoted to a similar problem 
- the problem of damping of the UHE neutrinos in the 
cosmic relic neutrino background, for example [2T ] [22 ] l23 l 



where a[s, Q^) is included in the integral over uj since 
it depends on the angle that also enters into the mo- 
mentum space element d'^o; (in other words, the parton 
originates from the photon and the momenta of these par- 
ticles are assumed to point in the same direction). The 
integration over e is performed over kinematically allowed 
parton energies. One can see that in pT]) . additionally to 
the familiar Planck's distribution 2/[e'^/^'^ — 1], we have 
introduced the function cpiejuj, Q'^)/lo to describe photon 
splitting into quark-antiquark pairs in a fashion analo- 
gous to the parton model. Actually, using that d'^oj — 
2TTLo^du) smOdO and making the change of variables ac- 
cording to (HI) (x = e/fcT, y = oj/kT) in ^ lead to 

7(i?)-^ j Ax{<j{xs,Q^)) 

X r^n(y)g-f^,g2V (12) 

Jx V \y J 

where the brackets (...) indicate the standard integra- 
tion over the polar angle 9, s — 2EkT{l—cos9), uq = N/V 
is the number density of the CMB photons (see 

Comparing ([T2]) with ([6]) we introduce the cross section 
averaged over the parton densities 

a{E) ^IJdx {a{xs, Q^)) q{x, Q^), (13) 

which allows to rewrite the damping rate in the ordi- 
nary form 

7(£;) = a{E)no. (14) 
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Fig. 4. Dependence of the survival probabihty of an UHE 
electron neutrino emitted at Zg — 300 (upper panel) and 500 
(lower panel) on its present day energy E. The dashed curve 
is obtained under the assumption that the only contributing 
channel is Ue'y W^X, while the solid one includes also the 
contribution of Ue'y ^ e~X 



For the reaction (El) one obtains 



dx ( 



(15) 



where cr,ye(s) is the cross section of the resonant sub- 
process z^eC^ W^. The integration over 6 includes the 
lepton distribution function since the latter, in this case, 
also depends on s (we deal with a Drell-Yan-like process 

Knowing ^{E), we can find the survival probability of 
an UHE neutrino traveling through the CMB radiation. 
For a constant damping, the probability as a function of 
E and the propagation time r is |23| 



P{E,t) 



(16) 



In fact, the travel without interaction can be so long 
that many cosmological quantities will have enough time 
to significantly change due to the expansion of the uni- 
verse (for example temperature of the CMB and therefore 
no). The neutrino energy will also get shifted to smaller 
values because of the expansion, just as in the case of light. 
This means that for a realistic situation the damping rate 
depends on time. In order to take this dependence into 
account it is convenient to use the redshift z which relates 
the present day values of the neutrino energy E and the 
CMB temperature T to those in the past. As we look back 
in time, these quantities increase with redshift as £'(l-f-z), 
T(l + z). We can then generalize (fT6|) by integrating the 



damping rate over all redshifts from now {z = 0) back to 
the neutrino source position Zg P^I^M^ : 



P{E, Zs) = exp 



H{z)il + z) 



^{E{l + z)) 



, (17) 



where dr = —dz/[H{z)(l + z)] with H{z) being the 
Hubble parameter (see e.g. [5^). Note that the CMB tem- 
perature entering into (jl7p (though not explicitly indi- 
cated) is understood also to scale as T{1 + z). We take 
H{z) = Ho^yO.S{^ + z)^ + 0.7 23J with Hq = 0.787 x 
10-28 cm-1 [27]. 

We first calculate the survival probability of the UHE 
neutrinos emitted at Zs — 10, 15 and 20 regarding the re- 
action ([8]) as the leading absorption channel. We parame- 
terize (Tye(s) in (HH) by the Breit-Wigner formula 



,(s) = 247r- 



r r 



p2 ' 



(18) 



where mw is the mass of the boson, F^e is the 
partial width for the decay W'^ — > I'sS^, and F is the 
total decay width of The corresponding results for 
the CMB temperature T = 2.725 K are shown in Fig. [3] 
(we set mw = 80.398 GeV, T^e = 0.230 GeV, F = 2.141 
GeV and a(m^) = 1/128 [27] in the positron distribu- 
tion function). One can see that the probability falls with 
increasing Zg (as intuitively expected) and the dip due to 
the resonance absorption becomes more distinct. 

We have also found the contribution of the reaction 
^ to the overall neutrino absorption. This turned out to 
be considerable at higher values of Zg. 

It has been stated above that in ([TS]) is fixed. 
Meanwhile we have to take all possible four-momentum 
transfers into account. To do this we have used the cross 



sections of the subprocesses of interest (fgd 
e~d) in the form 



dg2 



da„ 



2\2 



dQ2 



and 



(19) 



(20) 



where Gi? is the Fermi coupling constant. 

After replacing (t(x, Q^) in ([T^ by da^^/dQ^ and sub- 
stituting thus obtained damping rates into p7p we have 
integrated over (from 10 GeV^ to xs) besides the other 
integrations. The lower integration limit is chosen to be 10 
GeV^ since we are interested in neutrino-photon interac- 
tions which can produce final state hadrons. The corre- 
sponding results for T — 2.725 K, Zg = 300 and 500 are 
shown in Fig. H (we set Gp = 1.16637 x 10"^ GeV'^, 
and a = 1/137.035 in the quark distribution func- 
tions). One can see that the contribution of the channel 
([9]) grows with the neutrino energy E as well as with the 
distance from the source Zs ■ It is also notable that almost 
all the neutrinos emitted at z^ = 500 with energies about 
E{1 + Zs) = 25 X 10^4 eV would be absorbed by the CMB. 
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4 Conclusions 

We have discussed how parton distributions can be at- 
tributed to a photon gas. We have regarded the CMB ra- 
diation as a particular case and noted that its parton con- 
tent could be resolved by UHE cosmic neutrinos. We have 
considered the possibility of absorption of the UHE neu- 
trinos by the CMB via the channels ^ and Accord- 
ing to our calculations the universe turns out to be quite 
opaque to the neutrinos due to the presence of the CMB 
radiation. For example, more than 10 % of the UHE neu- 
trinos emitted at Zg = 15 with the initial energy 10^^ eV 
would be absorbed. Moreover, the absorption process is 
expected to be accompanied by hadron production thus 
providing a mechanism of generation of UHE cosmic rays. 
It is interesting to compare our results with the predictions 
of the survaval probability of the UHE neutrinos travers- 
ing the relic neutrino background (Ct^B) (see e.g. [2Tl[22l 
[^51l24j ). It is essential that the Ci^B has never been de- 
tected directly and there are only theoretical estimations 
of its parameters (such as temperature) while the proper- 
ties of the CMB have been well established experimentally. 
Therefore in studying the Ci^B by observing the resonant 
annihilation of the UHE neutrinos on the background an- 
tineutrinos through the reaction I'D one has to be 

able to distinguish between signals from this reaction and 
the interactions. This is possible, at least in principle. 
We have taken into account only the contribution of the 
charged current deep inelastic neutrino-photon scattering 
since in such a process, apart from the hadrons in the fi- 
nal state, the charged lepton is produced (in our case this 
is the electron). The latter being observed in correlation 
with the hadrons (or their decay products) may serve as 
an additional signature of the UHE neutrinos. 

I thank F. F. Karpeshin for highlighting a few important points 
on which I have focused more attention. This work was sup- 
ported in part by the Russian Foundation for Basic Research 
(grant 06-02-16135). 
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